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Abstract — In this paper, we propose a model for energy 
cooperation between cellular base stations (BSs) with individual 
renewable energy sources, limited energy storages and con- 
nected by resistive power lines for energy sharing. When the 
renewable energy profile and energy demand profile at all BSs 
are deterministic or known ahead of time, we show that the 
optimal energy cooperation policy for the BSs can be found 
by solving a linear program. We show the benefits of energy 
cooperation in this regime. When the renewable energy and 
demand profiles are stochastic and only causally known at 
the BSs, we propose an online energy cooperation algorithm 
and show the optimality properties of this algorithm under 
certain conditions. Furthermore, the energy-saving performances 
of the developed offline and online algorithms are compared by 
simulations, and the effect of the availability of energy state 
information (ESI) on the performance gains of the BSs' energy 
cooperation is investigated. 

I. Introduction 

In recent years, motivated by environmental concerns and 
energy cost saving considerations, telecommunication service 
providers have started considering the deployment of renew- 
able energy sources, such as solar panels and wind turbines, 
to supplement conventional power in powering base stations 
(BSs). In some places where the conventional power grid is 
still under-developed, the deployment of renewable energy 
sources is more attractive due to the significantly higher 
costs, as compared to a developed city, in powering BSs 
using conventional power sources. Examples where such a 
scenario occurs include the deployment of BSs with renewable 
energy sources by Ericsson in Africa [1] and Huawei in 
Bangladesh m. 

Although renewable energy sources are attractive for the 
above reasons, they also suffer from significantly higher 
variability as compared to conventional energy sources. As 
a result, even in BSs that deploy renewable energy sources, 
conventional energy sources, such as diesel generators or the 
power grid, is still required to compensate for the variability 
of the renewable energy sources. A method of mitigating 
the variability of renewable energy sources is through energy 
storage means such as fuel cells and batteries. Energy storage, 
however, is very costly to deploy and therefore, the amount 
of storage available at BSs will usually be quite limited. A 
key consideration in deploying BSs with renewable energy 
sources is minimization of the amount of conventional energy 
consumed, because it is only then cost-effective to deploy 
renewable energy sources and storage. A survey of issues in- 
volved in deploying renewable energy sources in BSs is given 



in Is). Related work on deploying renewable energy sources in 
smart grids, not necessarily constrained to a communications 
systems setup can be found in lH, ||5l and the references 
therein. 

In this paper, we consider mitigating the variability of 
renewable energy sources through geographical diversity. We 
consider the case when two or more BSs are connected by 
power lines so as to allow for transfer of energy between 
each other. A transfer of energy between two BSs allow 
for one that has excess of energy to compensate the other 
that has a deficit due to either higher demand of the users 
connected to the BS, or lower generation of renewable energy. 
We analyze the reduction in conventional energy needed to 
power the BSs if they are allowed to transfer energy, even 
when there is storage inefficiency and resistive power loss. 
We consider the availability of different information about 
the renewable energy sources and demand for our setting, 
and propose algorithms that take advantage of the energy 
cooperation between BSs and the information available to 
minimize their energy consumption from conventional sources. 

Another motivation for considering energy transfer comes 
from the possibility of using the power line as a backhaul 
link to enable coordinated multipoint transmission for cellular 
BSs 161. This results in an attractive dual use of the power line 
for both energy cooperation and communication cooperation. 

The rest of this paper is organized as follow. In Section |ll] 
we give formal definitions and description of our proposed 
energy cooperation model. In Section |III] we study the optimal 
offline energy cooperation policy for the case of deterministic 
renewable energy and demand profiles in which the future 
renewable energy and demand are known in advance. This 
setup, which has also been considered in energy harvesting 
based wireless communications Q, IS), models the scenario 
where we have good approximations of the renewable energy 
and demand profiles for the duration of interest and are 
willing to ignore small prediction errors. In Section |IV] we 
consider the general case of arbitrary renewable energy and 
demand profiles, and propose an online energy cooperation 
policy for this case. We analyze the optimality properties of 
this online policy under certain conditions, and compare its 
performance with the lower bound obtained by the offline 
policy via simulation. Finally, in Section [V] we conclude the 
paper and highlight some further directions in our extended 
version of this paper ||9l- 
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Fig. 1: System setup 

II. System Model 

In this paper, our focus will be on the case of two base 
stations, namely BS 1 and BS 2, with individual renewable 
energy generators, conventional energy sources, energy storage 
devices and connected with a power line. Our model, as 
depicted in Fig. [1] can be easily generalized to multiple (more 
than two) BSs, but we consider only the case of two BSs in 
this paper for simplicity. 

We consider a finite-horizon time-slotted system with slot 
index t, 1 < t < N, and N denoting the total number of slots 
under investigation. In the following, we define the elements 
of our energy cooperation model with two BSs, i.e. BS 1 and 
BS 2. We will use i G {1,2} to denote an element at the 
corresponding base station. 

A. Model Elements 

Renewable energy generated at BS i and time t: 

RE,{t) > 0. 
Demand at BS i and time t: DEi{t) > 0. 
Net energy generated at BS i and time t: Ei{t) = 

REi{t) — DEi{t). This quantity can be positive, representing 
a surplus, or negative, representing a deficit. 

Energy stored in BS i at time t: Si{t) > 0. To model 
limited storage constraint, we further assume Si{t) < 5max- 

Energy charged/discharged to/from storage at BS i and 
time t: c,{t) > /d,{t) > 0,d,(t) < s,(t). Intuitively, for 
given BS i and time t, there is at most one of Ci{t) and di{t) 
that is strictly positive, i.e. Ci{t) ■ di{t) = 0. 

Energy transfer from BS 1 (or 2) to BS 2 (or 1): xi2 (t) > 
(or X21 > 0). For a given time t, there is at most one of 
Xi2{t) and X2i{t) that is strictly positive, i.e. xi2{t) ■X2i{t) = 
0. 

Energy drawn from conventional energy source at BS i 
and time t: Wi{t) > 0. 

B. System Dynamics 

We require the following equations for storage dynamics to 

be satisfied: Si{t+1) — Si{t)+aci{t)~di{t). Here, < a < 1 



represents storage inefficiency, i.e. the energy lost in storage. 
As discussed earlier, we also require < Si (t) < S'max for all 
t. The combined storage dynamics and constraint lead us to 



the constraint: —Si{t) < aci{t) — di{t) < Sn 



s^{t). We 



also assume si(l) = S2(l) = 0. That is, there is no energy 
in storage at the initial time. Furthermore, the following two 
inequalities need to be satisfied at BS 1 and BS 2, respectively, 
in order to maintain their energy neutralization at each time 



t: 



Ei{t) + wi{t) - ci{t) + adi{t) - xi2{t) + Px2i{t) > 0, 



(1) 



E2{t) + W2it) - C2{t) + ad2{t) - X2l{t) + I3xi2{t) > 0. 

(2) 

Here, a again represents storage inefficiency and captures in 
this case, the inefficiency in drawing energy from storage. < 
/3 < 1 represents resistive loss in transferring energy from one 
BS to another ([T]i captures the constraint that any demand at 
time i at BS 1 has to be satisfied, by perhaps a combination of 
discharge from storage, transfer from BS 2 and conventional 
energy. Similarly, (|2) captures the energy balance requirement 
for BS 2. 

C. Control Policy and Objective Function 

In general, Ei{t) and i?2(i) can be modeled by 
a jointly distributed continuous stochastic process 
with a joint distribution F. Using vector notation, 
for any scalars yi{t),y2{t), . . . ,yn{t), we let 
y{t) = [yi{t),y2{t),---,yn{t)Y- Hence, we let 
s{t) ~ [si(t), S2(i)]"^ represent the state of our system 
at time t. Similarly, our control variables are the tuples 
{w{t) , c{t) , d{t) ^ xi2{t) , X2i{t)) . In general, these control 
variables at time t are functions of the past history, 
{{E{k)), I < k < t}, and the joint distribution F. A control 
policy TT is a sequence of these control variable^!]. That is, 

TT = {{wit),c{t),d{t),Xi2{t),X2lit)),l< t < N}. 

Then, the objective of our setup is to minimize the expected 
average conventional energy consumed. That is, we seek a 
control policy tt* that minimizes 
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where the expectation is taken with respect to the joint 
distribution F, and under the control policy tt* . 

Remark 2.1: Another valid cost criteria is to let iV — > cx) 
and minimize the long-run expected average conventional 
energy cost. That is, we wish to minimize 
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'The use of the symbol tt to represent a control policy is standard in the 
control/dynamic programming literature. With an abuse of notation, we will 
also be using the symbol tt to represent the number 3.14159... in our numerical 
simulations. It will be clear from context whether we are using the symbol tt 
for a control policy or the number. 



This criteria has the advantage of being insensitive to the 
starting state, but intuition about our model can be more 
easily obtained when N is finite. In this paper, we will restrict 
attention to finite N for simplicity. 

The optimal control policy for our model, as currently 
formulated, is open. In the rest of this paper, we will consider 
a number of special cases in which we can obtain some useful 
insight on this problem. 

III. Offline Algorithm with Deterministic Energy 
Profile 

The first restriction we make to this model is to consider a 
deterministic energy profile, with the net energy profile Ei (t) 
and E2{t) being known to both BSs for all t. In this case, our 
model reduces to the following linear program. 

Theorem 1: When the net energy profiles Ei{t) and E2{t) 
are deterministic and known to BS 1 and BS 2 for all t, the 
optimal control policy, tt*, is found by solving the following 
linear program. 

N 

miny^(ti;i(t) + W2{t)) 

TT ^ ^ 

t = l 

subject to (for 1 < t < iV) 

s{t + 1) s{t) + ac{t) - d{t), 

El {t) + wi (t) - ci (i) + adi (t) - xi2 (i) 

+ pX2lit) > 0, 

E2 (t) + W2 (t) - C2 (t) + ad2 [t] " X21 [t) 
+ pXi2{t) > 0, 

maxj 5 

s(l) = 0,s(2) = 0, 
c[t),d{t) > 0, 
d{t) < sit), 

Xl2{t),X2l{t) > 0. 

Proof: The reduction to the liner program follows from 
the assumption that the energy profiles are known for all t. 
In this case, the objective function simply reduces to the sum 
of the conventional energy required at each time t. Note that 
in the above problem, we do not explicitly put the constraints 
Ci{t)di{t) ~ Q, i — 1,2, and Xi2{t)x2i{t) ~ for any given 
t. However, it can be shown that the optimal solution of this 
problem always satisfies these constraints, and thus there is no 
loss of optimality in removing such constraints. ■ 

The assumption of deterministic energy profile models the 
case when the demands and renewable energy generated can 
be well approximated for 1 < t < N; that is, the case when the 
error in predicting the demand and renewable energy generated 
is small. Furthermore, it also allows us to gain insight into 
situations where it is beneficial for BSs to cooperate with 
each other. Intuitively, energy cooperation is helpful whenever 
the net energy generated at the two BSs are sufficiently 
uncorrected or anti-correlated, as will be shown next. 

To demonstrate the benefits of energy cooperation for two 
BSs, we model Ei{t) and E2{t) with the following energy 
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Fig. 2: Conventional energy consumed versus storage for 
different 6 

profiles. 

Ei{t) ^ Asm{ujt), (3) 
E2it) Asm{ujt + e). (4) 

Here, the correlation between the net energy profiles at BSs 
1 and 2 is measured by the phase shift 9. This approach 
of modeling correlation has been used in related context, 
such as in work on communications with energy harvesting 
devices |[T0|. 

Energy saving versus storage for different 6: We now show 
some simulation results on the energy saving versus storage for 
different values of 6. We set the following values: lu = 27r/24, 
A = 3,0<t< 239, 9 G {7r/4,7r/2,37r/4,7r}, /3 = 0.8 and 
a = 0.9. The results are plotted in Fig. |2] We compare the 
average cost {J2f'^Q{wi{t) + W2{t)) /2) against that of a single 
BS having the energy profile in ^ (plotted in green in the 
figure). 

As we can see from the figure, BSs' energy cooperation 
helps in general as the average cost per BS for the two 
cooperating BSs is lower than that of a single BS. As 9 varies 
from 7r/4 to tt, the cost per BS decreases since the energy 
profiles of the two BSs become more anti-correlated. 

As storage increases, it is also clear that the cost decreases, 
since more of the excess energy generated can be stored for 
later use, when there is a deficit. This storage benefit, however, 
decreases with 9 increasing to tt. Increasing 6' to tt signifies 
an increase in geographical diversity, resulting the ability to 
compensate deficit at one BS with excess from the other BS. 
When 9 — TT, there is little benefit from increasing storage. 

Energy savings versus 9 for fixed storage: To show the effect 
of 9 more clearly, we now keep the storage fixed at S'max = 
1 and vary 9 from to 27r. The rest of the parameters are 
kept fixed. In Fig. [3] we plot the percentage cost savings, 
relative to the energy cost of a single BS with the energy 
profile of (|3]l, against different values of 9. As we can see 
from the figure, the saved cost increases as 9 varies from to 
TT, at which point the energy profile of BS 2 is anti-correlated 




Fig. 3: Percentage energy saving versus 6 for S'max = 1 

with BS 1. This allows effective compensation through energy 
transfer between the two BSs. As 9 varies from tt to 2tt, the 
energy profile becomes highly correlated again, resulting in 
fewer opportunities to perform energy transfer between the 
two BSs. 

IV. Online Algorithm with Stochastic Energy 
Profile 

We now consider the more practical case when the net 
energy at both BSs are stochastic and not known ahead of time. 
We propose an online energy cooperation algorithm based 
on a greedy heuristic for minimizing conventional energy 
usage in Section IIV-AI We then analyze some properties of 
this algorithm in Section IIV-BI In particular, we state some 
optimality properties under specific energy profiles. Finally, 
in Section IIV-CI we provide simulation results on the perfor- 
mance comparison between the online algorithm versus the 
optimal offline algorithm proposed in Section |III] 

To describe the algorithm, we first assume a > and /3 > 
to avoid the complications of dealing with the case of no 
storage (a = 0) or no cooperation between BSs (/3 = 0). 

A. Greedy Online Algorithm 

Our greedy online algorithm for minimizing conventional 
energy consumed can be described by the actions BSs 1 and 
2 would take under the possible values of Ei{t) and E2{t) 
at time t as follows. Unless otherwise stated, we set all of 
w{t), c{t),d{t),x 12 (t) and X2i{t) equal to zeros for all t. 

Case 1: Ei{t) > and £'2(t) > 0. In this case, clearly, 
'w{t) = and both BSs try to store as much of the net energy 
as possible. That is, for i e {1,2}, we first carry out the 
following 

Ci{t) = min{(S'max - Si{t))/a,Ei{t)}, 

If both si{t) = S2{t) = S'max or si{t),S2{t) < 5max, tWs case 
terminates. Otherwise, if S2{t) < S'max and si{t) = S'max, BS 



1 transfers energy to BS 2 for storage. That is, we set 

Xi2{t) ^ Ei{t) - Ci{t), 

c'2{t) = min{/3xi2(t), (S^ax - S2(i))/"}, 
S2{t) ^ S2{t)+ac'^{t), 

C2{t)^C2(t)+c'2{t). 

Similarly, if if si(t) < Smax and S2{t) — Smax, the roles 
of BSs 1 and 2 in the above are reversed. 

Case 2: Ei{t) > and E2{t) < 0. In this case, BS 1 first 
transfers the net energy to BS 2 to compensate for the deficit. 
Hence, let 

xi2{t)^rRm{Ei{t),\E2{t)\/(3}, 

E'2{t)^E2{t)+l3xi2{t). 

Now, if i?2(i) = 0, we carry out the algorithm in Case 1 
with net energy profiles E[{t) = Ei{t) — xi2{t) and i?2(i)- 
If £'2(0 < 0, we compensate for the remaining deficit via 
storage at BS 2 first. We set 

d2{t)^Ta\xi{\E'^{t)\/a,S2{t)}, 
E'^{t)=E'2{t)+ad2{t). 

If £^2(0 = 0, this case is completed. Otherwise, we compen- 
sate from storage at BS 1. That is, we set 

=min{|i?^'(OI/M),Si(i)}, 
Xi2{t) Xi2{t) + adi{t), 
E'^'{t)=E'^{t)+aPd,it). 

Finally, if there is still a deficit remaining (£'2"(t) < 0), we 
compensate through conventional energy consumption and set 
W2{t) = \E'^'it)\. 

Case 3: Ei{t) < and E2{t) > 0. This case is symmetric 
to Case 2, with the roles of BSs 1 and 2 reversed. We therefore 
omit the description of the algorithm here. 

Case 4: Ei{t) < and E2{t) < 0. In this case, each BS 
compensates using individual storage first, before helping the 
other. That is, for i G {1, 2}, we let 

d^it) =mm{s^it),\E,{t)\/a}, 
E[{t)=E,{t)+ad,{t), 
Si{t) ^ s,{t) - dt{t). 

If either E[{t) > or E'^it) > 0, the algorithm reduces to 
the first three cases with net energy profiles being E'i{t) and 
E'^{t). If both E[{t) < and E'^{t) < 0, we compensate with 
conventional energy generation and set Wi{t) = 

For the case of no storage (a = 0), the greedy algorithm is 
modified in the obvious manner by not sending any energy to 
storage or drawing energy from storage. It is easy to show that 
this modified greedy algorithm is optimal for arbitrary energy 
profiles. 

For the case of no cooperation (/3 = 0), the greedy algorithm 
is again modified in the obvious manner by not requiring 
transfers between two BSs. Optimality of the greedy algorithm 
for this case will be discussed in the next subsection. 



B. Optimality Properties 

Although the greedy algorithm is a conceptually simple one, 
it has several optimality properties that we now analyze. Due 
to space constraints, we will only provide proof sketches of 
the propositions below, deferring the full proof to the extended 
version [[91 • 

Proposition 1: If /3 ~ or /3 = 1, the greedy algorithm is 
optimal for arbitrary energy profiles. 

Proof: When /3 = 1, the system reduces to a single 
BS with E{t) Ei{t) + E2{t) and < s{t) < 25max- 
The optimality of the greedy algorithm can then be inferred 
from Is] Theorem 1]. 

For the case of /3 = 0, no cooperation between the two BSs 
is possible. Optimality of the modified greedy algorithm for 
this case follows again from the fact that the greedy algorithm 
is optimal for individual BSs ||5] Theorem 1]. ■ 

We now proceed to analyze the greedy algorithm for the 
non-boundary cases of < /? < 1. It will be useful to define 
the following quantities. Define the unnormalized cost-to-go 
function under policy tt at time t and state s{t) as 



N 



J,(s(fc))=E^l^«;,(i) 



(5) 



^t=k 



We also denote the optimal cost-to-go function under the 
optimal policy tt* as J^r* {s{t)). For our setting, J^(s{t)) has a 
number of useful properties, which we now state. The results 
that follow in the rest of this section hold for any stochastic 
net energy profiles. Therefore, in our proofs, we will suppress 
the dependence of the control policy on the joint distribution 
of Ei{t) and E2it). 

Proposition 2: Suppose s'{k) > s{k) component- wise. 
Then, 



J,. (s(fc)) < J^, (s'(fc)) + (s'(fc) - s{k)) 



(6) 



Proof Sketch: The proof follows from the observation 
that a system with s(t) < s'{t) component-wise can mimic 
the optimal policy of a system at state s'{t) using conventional 
energy. 

Let 7r*(s'(fc)) denote an optimal policy when the state is 
at s'(fc), and {w* (t) , c* (t) , d* (t) , xl^it) , X2iit)) denote the 
control variables induced by the energy profile and optimal 
policy for t > k. We also use s*{t) to denote the evolution of 
the state under the optimal policy, starting from s'{k). 

Let 7r(s(/c)) denote a control policy when the state is at 
s{k), and {w{t),c{t),d{t),xi2{t),X2i{t)) denote the control 
variables induced by the energy profile and optimal policy for 
t > k. We use s{t) to denote the evolution of the state under 
the policy, starting from s{k). Now, we set tt = tt* except 
when d*{t) > Si{t) for any i £ {1,2}. In this case, we set 

d^it) = (7) 
w,{t)=w*{t)+a{d*{t)~s,{t)). (8) 

It is easy to show that the policy tt satisfies the constraints 
of the system and the maximum additional amount of con- 
ventional energy required, with respect to tt*, is at most 
al^(.s'(t) ■ 



The same proof strategy for a system at a different storage state 
s mimicking the optimal policy of the same system at storage 
state s' can be used to prove the following two propositions. 
We omit their proofs due to space constraint. 

Proposition 3: It is optimal to store excess energy at each 
of BS 1 and BS 2 first if there is still storage available, rather 
than to transfer the energy between them for storage. More 
concretely, suppose A > units of energy is available at BS 
1 at time t — k with si (fc) + a A < S'max and S2 (fc) + a/3 A < 
Sniax, then 



si(fc) + aA 

S2{k) 



si(fc) 
S2{k) + aPA 



Proposition 4: If /3 > a, then energy transfer is al- 
ways optimal. That is, if Ei{k) > > E2{k), sending 
min{|£;2(fc)|//3,i^i(fc)} units of energy from BS 1 to BS 2 at 
time t — k is an optimal policy. 

Remarks on Propositions ^ to ^ The above three propo- 
sitions show formally some intuitive aspects of energy co- 
operation that we would expect in such a system, and also 
certain optimality aspects of the greedy algorithm. It is obvious 
that a system starting with higher stored energy states has a 
lower optimum cost. Proposition |2] quantifies the maximum 
additional cost incurred by a system starting from a lower 
storage state. Proposition |3] shows the intuitively obvious fact 
that it is better to store energy locally rather than to store the 
energy at storage of the other BS. Proposition |4] formalizes the 
notion that if (3 > a, then it is more cost efficient to transfer 
energy to help the other BS rather than to store energy for 
future use, since the proportional loss in energy storage (1 — a) 
is higher than that in energy transmission (1 — /3). 

Proposition 5: If Ei (t) > for all 1 < t < iV and /3 > a, 
then the greedy policy is optimal. 

Proof Sketch: At each time t, there are two cases to 
consider, the case when E2 (t) > and the case when E2 (t) < 
0. In the first case, the greedy policy is optimal and follows 
straightforwardly from Proposition |3] 

In the second case, recall from Proposition |4] that energy 
transfer is an optimal strategy. If Ei{t) > \E2{t)\//3, then 
we reduce the problem back to the first case with E'i{t) = 
Ei{t) -\E2{t)\/l3 and E'2{t) = 0. For the case when there 
is still a deficit at BS 2 even after energy transfer (Ei{t) < 
\E2{t)\/f3), the optimality of the greedy algorithm is proved by 
first showing that it is optimal to minimize conventional energy 
usage at time t using arguments similar to that in Proposition |2] 
and the assumption that Ei (t) > 0, and then by showing that it 
is optimal to discharge storage at BS 2 first before discharging 
that at BS L ■ 

The condition that /3 > a can be relaxed, if more assump- 
tions can be made about the energy profile E2{t). 

Proposition 6: If Ei{t) > and E2{t) < for all 1 < t < 
N, then the greedy policy is optimal. 

Proof Sketch: The proof follows the same idea in proving 
the second case of Proposition |5] with the assumption that 
E2{t) is always non-positive instead of the condition (3 > a. 
When ^ < a, it can be more efficient in general to compensate 
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Fig. 4: Percentage energy loss of greedy online algorithm w.rt. 
offline algorithm versus 

for any deficit from local storage at BS 2 than to transfer 
energy from BS 1, since storage at BS 2 can be recharged 
more efficiently in future time steps with /3 < a. However, 
when E2{t) is always non-positive, any charging of storage at 
BS 2 has to come from BS 1, which incurs a proportional loss 
of 1 — aP >1 — a, resulting in any discharging or charging 
being less efficient than energy transfer. ■ 

C. Numerical Results 

We now compare the greedy online algorithm to the optimal 
offline algorithm proposed in Section Due to the lack of 
future information on the energy profiles, the greedy algorithm 
clearly cannot do as well as the offline algorithm, except 
under conditions discussed in Section IIV-BI We compare the 
performance differences between the two algorithms when the 
conditions in Section lTV-Bl are not satisfied. We adopt the same 
simulation setting as in Section Hill 

Fig. m shows the performance gap between the greedy 
online algorithm and the optimal offline algorithm that has 
access to the entire energy profile. Somewhat surprisingly, the 
greedy online algorithm suffers only a small loss (maximum of 
2.76%) compared to the offline algorithm under the sinusoidal 
energy profiles assumed in (|3} and (|4). Note that the vertical 
axis shows the loss of the greedy online algorithm, which is 
the percentage increase in energy consumption with respect 
to the optimal offline algorithm over the entire time horizon 
of N . Furthermore, when the energy profiles at the two BSs 
are anti-correlated {0 = tt), or highly positively correlated {6 
small or close to 27r), the percentage loss due to using an 
online algorithm is significantly lower, and the greedy online 
algorithm is almost as efficient as the offline algorithm. The 
offline algorithm, however, seems to be able to make more 
intelligence use of storage, resulting in generally higher energy 
saving when the the maximum amount of storage is increased. 
This suggests that greedy charge and discharge strategy may 
not be optimal, as can be expected. 



V. Conclusion and future directions 

In this paper, we have proposed a model for energy co- 
operation between two cellular BSs with individual renewable 
energy sources, limited storages, and a connecting power line. 
We consider two extreme cases. In the first case, we assume 
that the energy profile is known entirely for the duration of 
interest and have shown that the optimal policy can be found 
through solving a linear program. In the second scenario, 
we assume that we do not have any statistical information 
about the energy profile in future time steps. In this case, 
we have proposed a greedy online algorithm and analyzed 
the optimality properties of the greedy algorithm under some 
conditions. 

Instead of these two extremes, a more realistic scenario 
may be to assume that some, but not complete information, 
is known about the future energy profiles at the BSs. For 
example, it may be reasonable to assume that the energy 
profiles consist of a deterministic waveform in which small 
amount of random noise is added at each time step to model 
the prediction errors. In such a scenario, we can combine the 
proposed offline and online algorithms to arrive at a hybrid 
algorithm to leverage on the available information about the 
energy profile. Essentially, we can use the offline algorithm to 
determine the policy for the deterministic portion of the energy 
profile and then use the greedy algorithm to compensate for 
any differences induced by the random noise. We refer readers 
to our extended version [f9]| for work in this direction. 
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